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^ . The mass spectra of charmonia, bottomonia and Be mesons are calculated in the 

)^ \ framework of the QCD-motivated relativistic quark model based on the quasipoten- 

" I tial approach. The dynamics of heavy quarks and antiquarks is treated fully rela- 

' tivistically without application of the nonrelativistic l(? expansion. The known 

one-loop radiative corrections to the heavy quark potential are taken into account 
perturbatively. The heavy quarkonium masses are calculated up to rather high or- 
bital and radial excitations {L = 5, Ur = 5). On this basis the Regge trajectories are 
D I constructed both in the total angular momentum J and radial quantum number n^. 

' It is found that the daughter trajectories are almost linear and parallel, while parent 

trajectories exhibit some nonlinearity in the low mass region. Such nonlinearity is 
^ . most pronounced for bottomonia and is only marginal for charmonia. The obtained 

I results are compared with the available experimental data, and a possible interpre- 

tation of the new charmonium-like states above open charm production threshold is 



in 

O ' discussed. 

PACS numbers: 14.40.Pq, 12.39.Ki 

I. INTRODUCTION 

' At present a vast amount of experimental data on the heavy quarkonium spectroscopy 

^ has been accumulated pj] . The number of known states is constantly increasing. Thus, in 

the last eight years more than ten new charmonium-like states have been discovered (for 
a recent review see e.g. [3]). The total number of charmonium states, listed in the Parti- 
cle Data Group Listings [1], is 25 at present. Some of the new states (such as ric{2S), he, 
Xc2(2P), etc.) are the long-awaited ones, expected by quark models many years ago, while 
some others, with masses higher than the threshold of the open charm production, have 
narrow widths and unexpected decay properties There are theoretical indications that 
some of these new states could be the first manifestation of the existence of exotic hadrons 
(tetraquarks, molecules, hybrids etc.), which are expected to exist in QCD (see e.g. ^] and 
references therein). In order to explore such options, a comprehensive understanding of the 
heavy quarkonium spectroscopy up to rather high orbital and radial excitations is required. 
The LHCb Collaboration at the Large Hadron Collider (LHC) [3] plans to search for the 
bottom counterparts of the newly discovered charmonium-like states [5|. At present, the 
experimentally known bottomonium spectrum consists of 14 states [l[. Moreover, new ex- 
perimental data on the spectroscopy of Be mesons are expected from LHC . Therefore, the 
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investigation of the masses of the excited heavy quarkonia states represents an important and 
interesting problem. A reliable calculation of the masses of excited heavy quark- ant iquark 
states will allow to single out experimental candidates for the exotic multiquark states. 

In Ref. 0] we calculated the mass spectra of charmonia, bottomonia and Be mesons 
on the basis of a three-dimensional relativistic quasipotential wave equation with a QCD- 
motivated potential. For this calculation, a f^/c^ expansion up to the second order was 
used. Masses of several lowest orbital and radial excitations were obtained, mainly for 
the states lying under the open flavour production threshold. This investigation indicated 
that the charm quark is not heavy enough to be considered as nonrelativistic, especially for 
excited states Therefore, a reliable consideration of the highly excited charmonium states 
requires a completely relativistic treatment of the charmed quark without an expansion in 
its velocity. In this paper we extend the approach previously used for the investigations 
of light meson spectroscopy |8|, where relativistic quark dynamics was treated completely 
relativistically, to heavy quarkonia. Then the relativistic quasipotential, which determines 
the quark dynamics in heavy quarkonia, is an extremely non-local function in the coordinate 
space. In order to make it local, we replace the quark energies, entering the quark spinors, 
with the corresponding on-mass-shell energies. Such procedure makes the quasipotential 
local, but introduces a rather complicated nonlinear dependence on the bound state mass. 
The quasipotential equation with the complete relativistic potential can then be solved 
numerically using previously developed numerical methods In order to improve our 
description, leading radiative corrections to the heavy quark potential [Io| are also taken 
into account. Such corrections are suppressed by additional powers of a^, which are rather 
small for heavy quarkonia, and are known only in the framework of the f^/c^ expansion. 
Therefore we treat them perturbatively. The calculation of the masses of highly orbitally 
and radially excited states up to the fifths excitation is carried out. On this basis, the 
Regge trajectories for charmonia, bottomonia and Be mesons can be constructed both in 
the total angular momentum J and radial quantum number rir, and properties like linearity, 
parallelism and equidistance of these trajectories can be checked. There are reasons to expect 
that the parent Regge trajectories can be nonlinear due to the compactness of their ground 
and lowest excited states, which puts them in the region where both the linear confining and 
Coulomb parts of the quark-antiquark potential play an equally important role. Note also 
that the possibility of the assignment of the experimentally observed highly excited heavy 
quarkonium states to a particular Regge trajectory could help in determining their quantum 
numbers and elucidatine; their nature. 



In recent papers [ll|, |12| we investigated the possible interpretation of some new un- 
conventional charmonium-like states [2| as diquark-antidiquark tetraquarks. In particular, 
the relativistic dynamical calculation of the masses of such states was performed. Here we 
complement this study by calculating the spectrum of the highly excited conventional heavy 
quark-antiquark states in the same mass region. As a result, we will obtain a consistent pic- 
ture of the recently discovered heavy quarkonium states within the same relativistic quark 
model. 

A vast literature on the heavy quarkonium spectroscopy is now available. Therefore, we 
mostly refer to the recent comprehensive reviews j2|, [l3-15|, where the references to earlier 
review and original papers can be found. Recent investigations of highly excited heavy 



quarkonium states and their Regge trajectories can be found, e.g., in Refs. |16Hl8|. For very 



recent unquenched lattice QCD calculations see, e.g., Ref. [19[ and references therein 
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II. RELATIVISTIC QUARK MODEL 

In the relativistic quark model based on the quasipotential approach a meson is described 
by the wave function of the bound quark-antiquark state, which satisfies the quasipotential 
equation of the Schrodinger type [3] 
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\ 2jj,R 2^r) J (27r 

where the relativistic reduced mass is 



^A/(p) = / 7^^(p> q; M)vi/M(q), (1) 



and E'l, E2 are given by 



E-\ = , E-) = . (3) 

Here M = Ei + E2 is the meson mass, mi, 2 are the quark masses, and p is their relative 
momentum. In the center-of-mass system the relative momentum squared on mass shell 
reads 

^ ' 4M2 ■ ^ ^ 

The kernel V{p, q; M) in Eq. ([!]) is the quasipotential operator of the quark-antiquark in- 
teraction. It is constructed with the help of the off-mass-shell scattering amplitude, projected 
onto the positive energy states. Constructing the quasipotential of the quark-antiquark in- 
teraction, we have assumed that the effective interaction is the sum of the usual one-gluon 
exchange term with the mixture of long-range vector and scalar linear confining potentials, 
where the vector confining potential contains the Pauli interaction. The quasipotential is 
then defined by 

V{p, q; M) = Ui{p)u2{-p)V{p, q; M)u,{q)u2{-q), (5) 

with 

V(p, q; M) = ^asD^Ml^Y2 + V^,M^1^2;, + Konf (k), 

where is the QCD coupling constant, D^i, is the gluon propagator in the Coulomb gauge, 
7^ and u{p) are the Dirac matrices and spinors and k = p — q. 
The effective long-range vector vertex is given by 

r^(k)=7^ + ^a^,r, (6) 

where k is the Pauli interaction constant characterizing the anomalous chromomagnetic 
moment of quarks. Vector and scalar confining potentials in the nonrelativistic limit reduce 
to 

KonfW = il-e)iAr + B), 

Konf(0 = eiAr + B), (7) 

reproducing 

Konf (r) = Konf (0 + KLiir) = Ar + B, (8) 
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where e is the mixing coefficient. Therefore, in this hmit the Cornell-type potential is 
reproduced 

4 ry 

VMr) = -i,- + Ar + B. (9) 

All the model parameters have the same values as in our previous papers 0, 113] ■ The 
constituent quark masses ruu = md = 0.33 GeV, = 0.5 GeV, rric = 1.55 GeV, nib = 4.88 
GeV, and the parameters of the linear potential A = 0.18 GeV^ and B = —0.16 GeV have 
the usual values of quark models. The value of the mixing coefficient of vector and scalar 
confining potentials e = —1 has been determined from the consideration of charmonium ra- 
diative decays [3] and matching heavy quark effective theory (HQET). Finally, the universal 
Pauli interaction constant k = —1 has been fixed from the analysis of the fine splitting of 
heavy quarkonia ^Pj- states 0]. In this case, the long-range chromomagnetic interaction of 
quarks, which is proportional to {1 + k), vanishes in accordance with the fiux-tube model. 



III. QUARK-ANTIQUARK POTENTIAL 

The investigations of the heavy quark dynamics in heavy mesons indicate that the charm 
quark is not heavy enough to be considered as nonrelativistic. Indeed, estimates of the 
averaged velocity squared for the ground-state charmonium give the value (f^/c^) ~ 0.25. 
For excited charmonium states the (v'^/c^) values are even higher. Therefore, a reliable 
calculation of the charmonium spectroscopy requires a completely relativistic treatment of 
the charmed quark without an expansion in its velocity. 

The quasipotential ([5]) can in principal be used for arbitrary quark masses. The sub- 
stitution of the Dirac spinors into ([5]) results in an extremely nonlocal potential in the 
configuration space. Clearly, it is very hard to deal with such potentials without any ad- 
ditional approximations. In order to simplify the relativistic QQ potential, we make the 
following replacement in the Dirac spinors: 



ei,2(p) = V^ia + P' ->^i,2 (10) 

(see the discussion of this point in jo], [20[). This substitution makes the Fourier transforma- 
tion of the potential ([5]) local. 

The resulting QQ potential then reads 

V{r) = Vsi{r) + VsB{r), (11) 



where the explicit expression for the spin-independent Vsi(r) can be found in Ref. (8|. The 
structure of the spin-dependent potential is given by 



VsBir) = aL-S + b 



^(Si-r)(S2-r)-(Si-S2) 



- c Si • S2 + L • (Si - S2) + e (LSi) (LS2), 

(12) 

where L is the orbital angular momentum, Sj is the quark spin, S = Si -I-S2. The coefficients 
a, b, c, d and e are expressed through the corresponding derivatives of the Coulomb and 
confining potentials. Their explicit expressions are given in Ref. jsf. Since we also include 
the one-loop radiative corrections in our calculations, the strong coupling constant in the 
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static potential ([9]) should be replaced by the corrected constant ay [10 

,2 



) = ) 



1 + 



4 2 / TT 



(13) 



ai 



31 

y 
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-77, 



/3o 



11 - -n 



where nj is the number of flavours, /z is a renormalization scale and — 0.5772 is the Euler 
constant. 

The resulting quasipotential equation with the complete kernel (ITT!) is solved numerically 
without any approximations. The remaining one-loop radiative corrections, which are not 
included in the renormalized coupling constant (IT^ . are known only in the framework of 
the v/c expansion lol| ; therefore we treat them perturbatively. The additional contributions 
are the following [3, 
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(a) to the spin-independent part 

4/3oa2(/i2)ln(/ir) 



2tt 



ml 



4/3oag(/i^)ln(/xr) 
"3 27r r 



2mim2 

(b) to the spin-dependent part 



4 /3oag(/i^) 
3 



27r 



P 



,ln(/ir) (p ■ r)^ /ln(/ir) l' 



(14) 



6a 



Sb 
5c 



Sd 



+ 



rriQ / 3 vrr^ 



12 



+ 



mim2 3 vrr^ 
'1 1 \ 



/3o 



+ 7i? 



3 j + ^ In(yur) — 3 ^(^^77117712 r) 



/3o 



^-T77 + 7i?h7-o +ir ln(/xr) - - ln(Vmim2r) 



2 



2/',,2^ 



6 12 



27rr^ 



mi 



4a2(/i 



3mim2 vrr^ 
4 87ra2(/i2^ 



29 



'/3o 



/3o 



- t/3o + 7s — - 3 + — ln(/ir) - 3 ln(ym^m^r 



3mim2 37r 



12^° " y ~ 



mi — m2 ^ 1 mi + m2 



mi + m2 



mi — m2 



mi 



+ 



Stt 



'ln(/ir) + 7E\ , 21 



'ln(ymlm^r) + 7^;' 



r 

2/', ,2^1 



IGtt 



mt 



+ ( — + — 

,mi m2 



ml / 3 vrr^ 
27rr'^ 



7 /3o ^ 
3^12+^^ 



(15) 
(16) 



(17) 



y - 3 J + y ln(/ir) - 3 ln(^mim2 r) 



, "^2 
In — , 

mi 



(18) 



where for quantities quadratic in the momenta we use the Weyl ordering prescription [21 



{/(r)pyV = ^{{/(r),P^},F'}. 
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Since we consider only heavy quarks, for tlie dependence of the QCD couphng constant 
as(/i^) on the renormalization point /i^ we use the leading order result 

4:71 

"^^^'^ = /3oln(^VA^)- ^''^ 

In our numerical calculations we set the renormalization scale fi = 2mim2/{mi + and 
A = 0.169 GeV, which gives as = 0.315 for mi = m2 = rric (charmonium); as = 0.224 for 
mi = 7712 = 'nT'b (bottomonium) ; and = 0.286 for mi = m^, m2 = m;, {Be meson). ^ 
For the equal mass case (mi = m2 = m) the contribution of the annihilation diagrams 

^^'-^fi-ln2)m (20) 



3m2 \3 

which is of second order in as, must be added to the spin-spin interaction coefficient 6c in 
Eq. ([nD. 

Moreover, for the calculation of the bottomonium mass spectrum it is also necessary 
to take into account additional one-loop corrections due to the finite mass of the charm 
quark We considered these corrections within our model in Ref. 23 1 and found that 
they give contributions of a few MeV and weakly depend on the quantum numbers of the 
bottomonium states. The one-loop correction to the static QQ potential in QCD due to the 



finite c quark mass is given by [22, l23 



AVir, mj = --^^ [ln(v/^m,r) + 7^ + Ei{^77i,r)] , (21) 
y Tcr 

where 

^ , , dt , ~ 

Ei{x) = e — = -7s - Inx - 2^ 



t ^in-7il 



and Oo = 5.2. 



IV. RESULTS AND DISCUSSION 

A. Calculation of the masses 

We solve the quasipotential equation with the quasipotential (fTT]) . which nonperturba- 
tively accounts for the relativistic dynamics of both heavy quarks, numerically. Then we 
add the one- loop radiative corrections (fT^ - fl20l) and the additional one-loop correction for 
bottomonium due to the finite mass of the charmed quark (12T|) by using perturbation theory. 
The calculated masses of charmonia, bottomonia and the meson are given in Tables Hl lIIIl 
where n = n,, + 1, n,, is the radial quantum number, L, S and J are the quantum num- 
bers of the orbital, spin and total angular momenta, respectively. They are confronted with 
available experimental data from PDG [l|. Good agreement of our predictions and data is 
found. It is important to note that the nonperturbative relativistic treatment gives a better 



Note that these values of asip'^) were fixed in Ref. [f\ to optimize the fine and hyperfine splittings in 
charmonium and bottomonium (cf. Tables HI HI)) . 
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TABLE I: Charmonium mass spectrum (in MeV). 
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agreement with data than our previous heavy quarkonium mass spectrum calculation [7|], 
where only relativistic corrections up to j(? order were taken into account. However, the 
differences between former and new predictions are rather small for most of the low-lying 
states and become noticeable only for higher excitations, where relativistic effects turn out 
to be particularly important. 

The -Be meson states with J = L, given in Tables UTTt are mixtures of spin-triplet P-Z^l) 
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TABLE II: Bottomonium mass spectrum (in MeV). 
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and spin-singlet \^Ll) states: 

l^j) = |^Li)cos^+ |3LL)sin^, 

= -\^Ll) sine + \'^Ll) cos e, j = 1 = 1,2,?,... (22) 

where ^ is a mixing angle, and the primed state has the heavier mass. Such mixing occurs 
due to the nondiagonal spin-orbit and tensor terms in Eq. ( 1T2|) . The masses of physical 
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TABLE III: Be meson mass spectrum (in MeV). 
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states were obtained by diagonalizing the mixing matrix. The obtained values of the mixing 
angles 6 are close to the ones given in Ref. 

B. Regge trajectories 

In our analysis we calculated masses of both orbitally and radially excited heavy quarkonia 
up to rather high excitation numbers (L = 5 and rir = 5). This makes it possible to construct 
the Regge trajectories in the (J, M^) and {rir, M^) planes using the following definitions: 
a) (J, M^) Regge trajectory: 

J = aM^ + ao; (23) 
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FIG. 1: Parent and daughter {J,M'^) Regge trajectories for charmonium states with natural parity 
(P = (—1)''). Diamonds are predicted masses. Available experimental data are given by dots with 
particle names. The dashed line corresponds to a nonlinear fit for the parent trajectory. 



b) {rir, M^) Regge trajectory: 

rir = (3M^ + /3o, (24) 

where a, /3 are the slopes and aQ, (3q are the intercepts. The relations fl23|) and (|24|) arise in 
most models of quark confinement, but with difi'erent values of the slopes. 

In Figs. iniTO] we plot the Regge trajectories in the (J, M^) and (n^, M^) planes for char- 
monia, bottomonia and Be mesons. The masses calculated in our model are shown by 
diamonds. Available experimental data are given by dots with error bars and corresponding 
meson names. Straight lines were obtained by the fits of the calculated values. The 
fitted slopes and intercepts of the Regge trajectories are given in Tables [IV] and |Vl We see 
that the calculated charmonium masses fit nicely to the linear trajectories in both planes 
(maybe with the exception of the parent trajectories, where the J/ip and rjc mesons seem 
to have slightly lower masses). These trajectories are almost parallel and equidistant. In 
the bottomonium and Bf. meson sectors the situation is more complicated. The daughter 
trajectories, which involve both radially and orbitally excited states, turn out to be almost 
linear. On the other hand, the parent trajectories, which start from ground states, are ex- 
hibiting a nonlinear behaviour in the lower mass (excitation) region. Such nonlinearity is 
more pronounced in bottomonium. The origin of this nonlinearity can be easily understood, 
if one compares the mean radii of these states. The values of the mean square radii J{r'^) of 
charmonia. Be mesons and bottomonia, calculated in our model, are given in Table IVII The 
static potential of the quark-antiquark interaction is plotted in Fig. [11] (solid line). In this 
figure we also separately plot the contributions from linear confinement (dashed line) and 
of the modulus of the Coulomb potential (dotted line). As seen form Fig. [TTl the Coulomb 
potential dominates for distances less than 0.15 fm, while the confining potential is dominant 
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FIG. 2: Same as in Fig. [T]for charmonium states with unnatural parity iP = (— l)''"''^). 



for distances larger than 0.5 fm. In the intermediate region both potentials play an equally 
important role. Therefore the light mesons and charmonia (with the exception of the rjf. 
and J /ip which are in the intermediate region) have characteristic sizes which belong to the 
region, where the confining potential dominates in the interquark potential ([9]). This leads 
to the emergence of the linear Regge trajectories. Contrary, the ground and few first excited 
states of bottomonia and B^. mesons have smaller sizes and fall into the region, where the 
Coulomb part of the potential ([9]) gives an important contribution. As a result, the parent 
Regge trajectories of bottomonia and Bf. mesons are nonlinear, while the daughter trajec- 
tories (which fall into the region, where the confining potential is dominant) are still linear 
ones. In Refs. [l6|, 24] an interpolating formula between the limiting cases of pure Coulomb 
and linear interactions was proposed. It can be written as follows: 
(a) for the parent trajectory in the (J, M^) plane 



M^=\J 



(7T^ + ^°) 



(25) 



(b) for the J = 1 trajectory in the (n^, M^) plane 



rir — 



+ To /r. 



(26) 



where the parameters 7, r, 70, tq and 71, ri determine the slopes, intercepts and nonlinearity 
of the Regge trajectories, respectively. Their fitted values are given in Table IVIII The 
corresponding Regge trajectories are plotted in Figs. [TlfTUlby dashed lines. It is found that 
these nonlinear trajectories have the same slope 7 for the given quarkonium family, which 
is generally in agreement (but slightly higher) with the linear trajectory slopes a for the 
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FIG. 3: The (n^-jM^) Regge trajectories for vector (5-wave), tensor and vector (D-wave) charmo- 
nium states (from bottom to top). Notations are the same as in Fig. [TJ 



respective daughter trajectories given in Table IIVI We see that the nonUnearity of the 
charmonium Regge trajectories is almost negligible, and its account does not significantly 
improve the quality of the fit compared to the linear one. 

We can compare the slopes of linear Regge trajectories for heavy quarkonia obtained 
in this paper with our previous results for the slopes of Regge trajectories of light [8l] and 



heavy-light [25|] mesons. Such comparison shows that the slopes decrease rather fast with 
the growth of the quark masses: the slope a decreases from about 1.1 GeV~^ for light 
mesons, composed from m, d quarks and antiquarks, to about 0.24 GeV~^ for bottomonium. 
However, the difference between slopes of heavy-light [Qq) mesons and of heavy quarkonia 
{QQ') is not so dramatic. In fact, comparing the present Tables HVl IVl and the corresponding 
Tables 4, 5 of Ref. |25|], we see that the slopes for charmonia and D, Dg mesons have very 
close values, and the same is true also for bottomonia, and B, Bg mesons. This might 
indicate that the slope of the meson Regge trajectory is mainly determined by the mass of 
the heaviest quark mq. The dependence of the Regge slopes a and /3 on trq has in both 
planes with rather good accuracy the same simple form: oc 1/ ^JnTq. 

From the comparison of the slopes in Tables llVt IVl we see that the a values are system- 
atically larger than the fi ones. The ratio of their mean values is about 1.4 for charmonia, 
bottomonia and B^ mesons. This value of the ratio is the same as for the charmed and 
bottom mesons 25|, but slightly larger than for light mesons for which a/ (3 was found 
to be about 1.3. 



C. Comparison with experiment 

In Tables [Il lIIII we compare our predictions for the heavy quarkonium masses with the 
available experimental data [l|. We find that all states below the open flavour production 
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FIG. 4: Same as in Fig. [T]for bottomonium states with natural parity. 



thresholds are well described by our model, the difference between predicted and measured 
masses does not exceed a few MeV. For higher excited states, which are above this thresh- 
old, most of the well-established conventional states (believed to be quark-antiquark ones) 
are also well described by our model, the difference between theory and experiment being 
somewhat larger, but it still does not exceed 20 MeV. ^ As it is seen from Figs. [TlfTO| these 
states fit to the corresponding Regge trajectories both in the (J, M^) and (nr.,M^) planes. 

We first discuss the recently found quarkonium states below the open flavour production 
threshold. The observation and measurement of the mass of the pseudoscalar ground state 
Vb 13] provides a significant information about the spin-spin interaction in heavy quarkonia. 
The averaged bottomonium hyperfine splitting measured in T{3S) — )■ nJlS)^ and T{2S) 
r7fe(l^)7 decays is AMMs(r/b) = Mr{is) - = 69.3 ± 2.8 MeV [H, 13. Very recently the 

Belle Collaboration [29| reported the first observation of the radiative transition hh{lP) 
r]b{lS)'y. The measured rjbilS) mass is 9401.0 ± 1.9l2!4 MeV and the hyperfine splitting 
AMMsiVb) = 59.3±1.9;i;4 MeV Q. Our prediction for this splitting, AMMsiVb) = 62 MeV, 
is in good agreement with the experimental values. ^ Note that our model correctly predicts 
the branching ratios of the corresponding radiative decays 0] ■ For the better understanding 
of the hyperfine interaction in heavy mesons it will be very interesting to measure the 
mass of the vector i?*-meson, which consists of two heavy quarks of different flavours, and 
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Note that hadron loop effects [26|, |27[ can lead to mass shifts and state mixings. As shown in Ref. [27[, 
the loop mass shifts can be absorbed by a change of the valence quark model parameters. However, the 
opening of the new threshold could provide a larger mass shifts to the nearby quarkonium states. 
Almost the same value of AMi^[s{t]ij) = 60 MeV was predicted by us in Ref. [7|, while most of other 
theoretical predictions, especially the ones based on perturbative calculations, gave significantly lower 



central values, e.g., 41 ± 14 MeV [30 1. 
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FIG. 5: Same as in Fig. [T]for bottomonium states with unnatural parity. 

determine the related hyperfine sphtting. Our model correctly predicts the pseudoscalar 
Bf. meson mass and gives for the hyperfine splitting the value of AMhf(-Bc) = 61 MeV. 
Note that the LHCb Collaboration very recently measured the B^. meson mass, and their 
preliminary result is M(_B+) = 6268.0 ±4.0(stat) ± 0.6(syst) MeV, leading to the improved 
average M(S+)^^p = 6272.95 ± 5.17 MeV [gj. 

Another important experimental test of the structure of the spin splittings in heavy 
quarkonia comes from the measurement of the masses of the spin-singlet P-levels first in 
charmonium hc{lP) jUj and very recently in bottomonium hb{lP) and hb{2P) j32[. The 
measured masses of these states almost coincide with the spin-averaged centroid of the 
triplet states {M{^Pj)) = [M(xqo) + 3M(xqi) +5M(xq2)]/9. The hyperfine mass splittings 
AMhfe(raP) = {M{n^Pj))-M{n^Pi) in bottomonium are found to be AMhfs(lP) = (1.62 ± 
1.52) MeV and AMhfs(2P) = (0.48^11^) MeV This observation indicates that the 

spin-spin contribution is negligible for P-levels, and thus shows the vanishing of the long- 
range chromomagnetic interaction in heavy quarkonia. In our model this is the result of the 
choice of the value of the long-range chromomagnetic quark moment k = — 1. Note that our 
original predictions [Tj] for the spin-singlet masses are confirmed by these measurements. 

The recently observed T{1^D2) state is the only D-wave state found below the threshold 
of open flavour production. Our prediction for its mass is consistent with the measured 
value. It will be interesting to observe other T{1D) states in order to test further our 
understanding of spin-orbit and spin-spin interactions in heavy quarkonia. 

Next we discuss the observed states above the open flavour production threshold. The 
most well-established states are the vector 1 states. For charmonium PDG [l| lists seven 
such states: ^(3770), ^(4040), V^(4160), X(4260), X(4360), 7/;(4415) and X(4660), from 
which only the ip states are included in the PDG Summary Tables [ij. These states are 
believed to be ordinary cc charmonium (with isospin / = 0). They are well described 
by our model (see Table [I]): '?/'(4040) and '?/'(4415) are the 3^ Si and 4^5*1 states, while 
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FIG. 6: Nonlinear Regge trajectories in the {J,M^) starting from vector, pseudoscalar and scalar 
bottomonium states (from bottom to top). 
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FIG. 7: Same as in Fig. [3] for bottomonium. 



'?/'(3770) and '?/'(4160) are the I'^Di and 'l^D^ states, respectively. These states fit well 
to the corresponding Regge trajectories (see Fig. [3]). On the other hand, the three new 
vector states X are considered as unexpected exotic states (their isospin is not determined 
experimentally). Indeed, we do not have any cc candidates for these states in Table [B 
Contrary, in Ref. [13] we have found that these states can be described in our model as 
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FIG. 9: Same as in Fig. [3] for the Be meson. 



tetraquarks composed from a diquark and antidiquark ([eg] [eg], q = u,d). In particular, the 
X(4260) and X(4660) states can be interpreted as the 1 states of such tetraquarks with 
a scalar diquark [cg]s=o ^^nd scalar antidiquark [eg] 5=0 in the relative IP- and 2P-states and 
predicted masses 4244 MeV and 4666 MeV, respectively [l^]. The X(4360) can be viewed 
as the 1 tetraquark with the axial vector diquark [cg]s=i and axial vector antidiquark 
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FIG. 10: Same as in Fig. [6] for the i?c meson. 
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FIG. 11: Static potential of the quark-antiquark interaction ([9]) without the constant term (sohd 
Une). Dashed line shows the hnear confining potential contribution, while dotted line corresponds 
to the modulus of the Coulomb potential. 



=1 in the relative IP-state, which mass is predicted to be 4350 MeV ^ 



Note that in Ref. 33| a different prescription for these states is used: X(4260), X(4360), ^'(4415) and 
X(4660) are assigned to 45", 31?, hS and 6S' charmonium cc states. This leads the authors to the conclusion 
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TABLE IV: Fitted parameters of the {J,M'^) parent and daughter Regge trajectories for heavy 



quarkonia and Be mesons with natural and unnatural parity. 



Trajectory 


natural parity 


unnatural parity 


a (GeV-2) 




a (GeV-2) 




cc 


T/ib 




'/c 




parent 






n 41 (J _|_ n noi 


_Q on _|_ n Q1 


first daughter 






n 470 _|_ Q ni cj 


_« _|_ n 24 


second daughter 


0.431 ± 0.036 


-6.08 ± 0.68 


0.414 ± 0.050 


-6.66 ± 0.92 


cc 


ACU 




ACi 




parent 


431 ± 016 


—5 07 ± 25 


461 ± 008 


—4 66 ± 12 


daughter 




1 .'ii _1_ W.fJO 




_tr 00 _|_ n 1 1 


bb 


T 




'lb 




parent 


212 ± 022 


— 1 8 5 ± 2 3 


184 + 024 


— 16.7 ± 2.5 


first daughter 


n 24fi + n 014 




n OQ4 _|_ n ni « 


_2T 5 + 17 


second daughter 


\J , ^ _1_ VJ.VJXW 


^ 1 . ± _1_ ± . ± 


n 040 1 n n-|4 




third daughter 


0.246 ± 0.027 


-26.6 ± 3.1 


0.241 ± 0.026 


-27 ± 3 


bb 


A/Ln 
AOO 








parent 


0.228 ± 0.021 


-22.3 ±2.2 


0.239 ±0.018 


-22.5 ± 1.9 


daughter 


0.254 ± 0.009 


-26.7 ± 1.0 


0.267 ± 0.006 


-27.1 ±0.7 


be 






Be 




parent 


0.254 ±0.018 


-9.38 ±0.88 


0.242 ±0.019 


-9.75 ±0.95 


first daughter 


0.291 ± 0.008 


-12.8 ±0.4 


0.278 ± 0.009 


-13.1 ±0.5 


second daughter 


0.270 ±0.010 


-13.3 ±0.4 


0.259 ±0.007 


-13.5 ±0.4 


be 






Bel 




parent 


0.265 ±0.013 


-12.0 ±0.6 


0.285 ± 0.007 


-13.0 ±0.4 


daughter 


0.275 ±0.014 


-13.9 ±0.8 


0.298 ± 0.008 


-15.3 ±0.4 



The three vector bottomonium states, T(10580), T(10860) and T(11020), observed above 
open bottom threshold [l[, are rather well described in our model as 4^5*1, 5^ Si and 6^ Si 
states (see Table [IT|) . The mass of T(11020) being somewhat higher than the experimental 
value. They fit to the corresponding Regge trajectory in Fig. [71 

The only experimentally established 2P charmonium state is Xc2{'2P) which mass is pre- 
dicted slightly higher (by about 20 MeV) in our model. From Tabled we see that the exotic 
state X(3872) cannot be described as the l"*""*" 2^Pi cc state or the 2~+ 1^D2 cc state. If this 
state belonged to either 2P or ID multiplets, this could signal a large fine splitting in these 
multiplets, since the X(3872) mass is 55 MeV below Xc2(2P) and 100 MeV above V'CSTTO]^ As 
we see from Table [U our model does not support such large fine splittings. In Ref. [111. Il2 



we 



argued that X(3872) can be considered as the 1++ ground state tetraquark, composed from 
the scalar and axial vector diquark and antidi quark (([cg]5=o[c^]5=i ± [cq]s=i[cq]s=o) / V^)), 
which mass is predicted to be 3871 MeV. 



that the quark interaction potential should be screened at large distances. 
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TABLE V: Fitted parameters of the {nr,M'^) Regge trajectories for heavy quarkonia and Be 



mesons. 


Meson 


P (GeV-2) 


Po 


Meson 


P (GeV-2) 


Po 


cc 






cc 






Vc 


0.287 ±0.011 


-2.62 ±0.18 




0.297 ±0.010 


-2.89 ±0.16 


XcO 


0.288 ± 0.003 


-3.34 ±0.06 


Xcl 


0.301 ±0.011 


-3.67 ±0.19 


Xc2 


0.301 ±0.011 


-3.76 ±0.019 


he 


0.298 ±0.010 


-3.68 ±0.17 




0.325 ± 0.006 


-4.62 ±0.11 




0.315 ±0.003 


-4.53 ±0.06 




0.311 ± 0.002 


-4.53 ± 0.04 




0.317 ±0.003 


-4.53 ±0.06 


bb 






bb 






Vb 


0.151 ± 0.013 


-13.7 ± 1.4 


T 


0.153 ± 0.012 


-14.1 ±1.3 


XbO 


0.158 ± 0.008 


-15.4 ±0.9 


Xbl 


0.159 ± 0.007 


-15.7 ±0.8 


Xb2 


0.161 ±0.007 


-15.9 ±0.7 


hb 


0.161 ±0.007 


-15.8 ±0.8 




0.178 ± 0.002 


-18.4 ±0.3 




0.178 ± 0.002 


-18.4 ±0.3 




0.178 ± 0.003 


-18.4 ±0.3 




0.178 ± 0.002 


-18.4 ±0.3 


be 






be 






Be 


0.172 ± 0.008 


-6.88 ± 0.39 


B*e 


0.175 ± 0.008 


-7.15 ±0.39 


Bco 


0.184 ±0.001 


-8.28 ±0.07 


Be2 


0.185 ± 0.001 


-8.48 ± 0.07 




0.190 ±0.002 


-9. 10 ±0.12 




0.188 ±0.002 


-9.27 ± 0.11 



TABLE VI: Mean square radii ^ {r'^) for the spin-singlet ground and excited states of charmonia, 
Be mesons and bottomonia (in fm). 



State 








IS 


0.37 


0.33 


0.22 


IP 


0.59 


0.53 


0.41 


2S 


0.71 


0.63 


0.50 


ID 


0.74 


0.67 


0.54 


2P 


0.87 


0.79 


0.65 


IF 


0.87 


0.79 


0.65 


3S 


0.94 


0.87 


0.72 


IG 


0.98 


0.89 


0.75 


2D 


0.99 


0.90 


0.76 


IH 


1.08 


0.99 


0.85 


3P 


1.09 


0.99 


0.84 


2F 


1.09 


0.99 


0.85 


4S 


1.16 


1.05 


0.90 


3D 


1.18 


1.08 


0.94 


4:P 


1.26 


1.16 


1.01 


5S 


1.32 


1.21 


1.07 


6S 


1.46 




1.22 
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TABLE VII: Fitted parameters of the nonlinear Regge trajectories for heavy quarkonia and Be 
mesons. 



Meson 


7 (GeV-2) 


70 


71 


T (GeV-2) 


To 


n 


T 


0.33 


32.2 


32.3 


0.22 


22.2 


10.1 


% 


0.33 


32.9 


15.0 








X60 


0.33 


33.7 


6.57 








B* 


0.32 


13.3 


12.5 


0.21 


9.25 


4.00 


Be 


0.32 


14.2 


6.21 








Bco 


0.32 


15.1 


2.99 








J/ij 


0.48 


4.25 


5.47 


0.31 


3.19 


0.82 




0.48 


5.19 


3.56 









As we see from Table [U the X(4160) and X(4350) can be attributed from the point of 
view of the mass value and charge parity C = + to the pseudo tensor 2~+ spin-singlet 
2^D2 and tensor 2~^^ spin-triplet 3^P2 charmonium states, respectively. They fit well to the 
corresponding Regge trajectories in Figs. [TH3l 



The X(4140) state, observed by CDF in — )■ K^(j)J/ip decays [3J], can correspond 
in our model to the scalar 0++ charmed-strange diquark-antidiquark [cs]s'=i[cs]5=i ground 
state, which predicted mass is 4110 MeV, or the axial vector l"*""*" one ([cs]5=o[cs]s'=i + 
[cs]s=i[cs]s=o)/V2) with calculated mass 4113 MeV [III, [III. 



Two of the three charmonium-like charged X states reported by Belle [35|, which are 
explicitly exotic, can be interpreted in our model as tetraquark states. We do not have 
tetraquark candidates for the X(4040)^ structure, while the X(4250)^ can be consid- 
ered as the charged partner of the 1~ IP state [cm] 5=0 [erf] 5=0 or as the 0~ IP state of 
the ([cu]5'=o[crf]s'=i + [cu\^i[cd]s=o) / \/2) tetraquark with predicted masses 4244 MeV and 
4267 MeV, respectively [l2|. The X(4430)+ could be the first radial (25) excitation of the 
1+ X(3872) tetraquark or the 0"*" 2S [cu]s=i[cd\s=i tetraquark, which have very close masses 
4431 MeV and 4434 MeV [11. 



Very recently the Belle Collaboration [36l] reported the observation of two charged 
bottomonium-like resonances, the Zfc(10610) and ^^(lOGSO), in the tt^T and n^hfy mass 
spectra close to the open-bottom {BB* and B*B*) production thresholds. The analysis of 
the charged pion angular distributions favours a = 1~*" spin-parity assignment for both 



states [36| . In this mass region we do not have any bottom diquark-antidiquark tetraquarks 
with such quantum numbers 13] . The possible interpretations of these exotic bottomonium- 
like states are discussed in Ref. 15 



As we see, a consistent picture of the excited quarkonium states emerges in our model. 
All well-established states and most of the states, which need additional experimental confir- 
mation, can be interpreted as excited quarkonium or diquark-antidiquark tetraquark states. 



V. CONCLUSIONS 

The mass spectra of charmonia, bottomonia and Be mesons were calculated in the frame- 
work of the QCD-motivated relativistic quark model based on the quasipotential approach. 
Highly radially and orbitally excited quarkonium states were considered. To achieve this 
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goal, we treated the dynamics of heavy quarks in quarkonia completely relativistically with- 
out application of the nonrelativistic w^/c^ expansion. The known one-loop radiative cor- 
rections were also taken into account in order to improve the agreement with experiment. 
The comparison of new results with the previous consideration within the v'^/c^ expansion 
[3] indicates that relativistic effects become significant with the increase of the excitation 
and are particularly important for charmonium. 

On this basis, the Regge trajectories of heavy quarkonia were constructed both in the 
(J, M^) and (ra^, M^) planes. A different behaviour of these trajectories was observed for 
parent and daughter trajectories. All daughter trajectories turn out to be almost linear 
and parallel, while parent trajectories exhibit some nonlinearity. Such nonlinearity occurs 
only in the vicinity of ground states and few lowest excitations and is mostly pronounced 
for bottomonia. For charmonia this nonlinearity is only marginal, and its account does not 
significantly improve the fit. It was shown that the masses of the excited states of heavy 
quarkonia are determined by the average distances between quarks larger than 0.5 fm, where 
the linear confining part of the quark- ant iquark interaction dominates. This leads to the 
emergence of almost linear Regge trajectories. On the other hand, a few lowest states 
have average sizes smaller than 0.5 fm and fall in the region, where both the Coulomb and 
confining potentials play an important role. As a result, the parent Regge trajectories exhibit 
a certain nonlinearity in this region. The parameters (slopes, intercepts and nonlinearity) 
of both linear and nonlinear Regge trajectories were determined. They were compared to 
the slopes of the hnear Regge trajectories of light [Sj and heavy-light |25j mesons calculated 
previously. It was found that the slope of the meson Regge trajectory is mainly determined 
by the mass of the heaviest quark mq. 

A detailed comparison of the calculated heavy quarkonium masses with available exper- 
imental data was carried out. It was found that all data for the states below open flavour 
production threshold are well reproduced in our model: the difference between predicted and 
measured masses does not exceed a few MeV. For higher excited states, which are above 
this threshold, most of the well-established conventional states are also well described by 
our approach, the difference between theory and experiment being somewhat larger, but still 
within 20 MeV. It was shown that these states fit well to the corresponding Regge trajec- 
tories. Other states, which have unexpected properties and are therefore believed to have 



an exotic origin, were also discussed. As it was shown in our previous calculation [11|, [12 
most of these states can be described as diquark-antidiquark tetraquarks. Therefore we have 
a self-consistent picture of the heavy quarkonium spectra. Future experimental studies of 
the yet unobserved conventional quarkonium states and a clarification of the nature and 
quantum numbers of the exotic quarkonium-like states will provide a further test of our 
model. 
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